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OO 

■ 1 Introduction 

O 

■ Routh's procedure, in its original form (as described in his treatise [13]), was concerned with 
eliminating from a Lagrangian problem the generalized velocities corresponding to so-called 
ignorable or cyclic coordinates. Let L be a Lagrangian on R n that does not explicitly depend 
on m of its base variables, say the coordinates 9 a . From the Euler-Lagrange equations for these 

K> , coordinates, 

' A ( —\ dL - n 

dt \W) ~d¥ ~~ ' 

we can immediately conclude that the functions dL/d6 a are constants, say 

dL 

these equations express the conservation of generalized momentum. Routh's idea is to solve 
these equations for the variables 9 a and to introduce what he calls the 'modified Lagrangian 
function', the restriction of the function 

L> = L-^6° 

de a 

to the level set where the momentum is 7r a . One can easily verify that the {n—m) Euler-Lagrange 
equations for the remaining variables x % can be rewritten as 

d (9L\ dL _ Q d (9L'\ dL' _ 

dt V dx l J dx l dt V dx l J dx % 
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For example, if the Lagrangian takes the form 



L(x,e,x,e) = lk ij {x)i i x j + k ia {x)x l e a + \k ab { x )e a e b - v(x), 

the conservation of momentum equations read ki a x % + k ab 8 = 7T a , and they can be solved for 
the variables 9 a if (k ab ) is a non-singular matrix. The modified Lagrangian function is 

L (x,x) = \{kij — k ki a kjf))x i? + k ki a ix b x — (V + ^k vr a 7r^), 

where k ab denotes a component of the matrix inverse to (k ab ), in the usual way. Clearly, the 
advantage of this technique is that the reduced equations in L' involve only the unknowns x 1 
and x l ; they can in principle be directly solved for the x l , and the 6 a may then be found (if 
required) from the momentum equation. 

A modern geometric interpretation of this reduction procedure can be found in e.g. [8] . The above 
Lagrangian L is of the form T — V, where the kinetic energy part is derived from a Riemannian 
metric (i.e. we are dealing with a so-called simple mechanical system). The function L is defined 
on the tangent manifold of a manifold of the form M = S x G (in this case R n ) and it is invariant 
under an Abelian Lie group G (in this case the group of translations R m ). The main feature 
of the procedure is that the modified Lagrangian function and its equations can be defined in 
terms of the coordinates on S only. However, to give the definition of the modified function an 
intrinsic meaning, we should define this function, from now on called the Routhian, rather as 
the restriction to a level set of momentum of 

FIT 

n = L-^(e a + Afx i ), 

with Af = k ab k ib , i.e. 

K(x,x) = - k ab k ia k jb )x i x i - (V + lk ab w a TT b ). 

The coefficients A" form a connection on the trivial principal bundle M = S x G — > S , usually 
called the mechanical connection, and 9 a + Afx 1 is in fact the vertical projection of the vector 
(x l ,9 a ). The (n — m) Euler-Lagrange equations in x % then become 

dtKdxO o,- 



where in the term on the right-hand side 



r>a 



dAf 3A« 



dxi dx l 

has a coordinate-free interpretation as the curvature of the connection. 

In [HKIO], Marsden et al. extended the above procedure to the case of simple mechanical systems 
with a non-Abelian symmetry group G and where the base manifold has a principal bundle 
structure M — * M/G. The procedure has recently been further extended to cover Lagrangian 
systems in general by Castrillon-Lopez pp. 

The most important contribution of our paper lies in the geometric formalism we will adopt. 
The bulk of the literature dealing with different types of reduction of Lagrangian systems has 
relied heavily on methods coming from the calculus of variations. In fact, as in e.g. [U [5], [TO] , the 
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reduced equations of motion are usually obtained by considering some reduced version of Hamil- 
ton's principle. Our method is different from those of other authors in that it doesn't involve 
consideration of variations. It is distinctively Lagrangian (as opposed to Hamiltonian) , and is 
based on the geometrical analysis of regular Lagrangian systems, where solutions of the Euler- 
Lagrange equations are interpreted as integral curves of an associated second-order differential 
equation field on the velocity phase space, that is, the tangent manifold of the configuration 
space. Consequently our derivation of Routh's equations is relatively straightforward and is 
a natural extension of that used by Routh in the classical case. In particular, we will show 
how Routh's equations can be derived directly from the Euler-Lagrange equations by choos- 
ing a suitable adapted frame, or equivalently by employing well-chosen quasi-velocities. This 
line of thinking has already provided some new insights into e.g. the geometry of second-order 
differential systems with symmetry [3]. 

We deal from the beginning with arbitrary Lagrangians, i.e. Lagrangians not necessarily of the 
form T-V. 

As in [TO] , we explain how solutions of the Euler-Lagrange equations with a fixed momentum can 
be reconstructed from solutions of the reduced equations. The method relies on the availability 
of a principal connection on an appropriate principal fibre bundle. We will introduce in fact two 
connections that serve the same purpose. 

We describe the basic features of our approach in Section 2. The reduction of a Lagrangian 
system to a level set of momentum is discussed in Section 3, and our generalization of Routh's 
procedure is explained there. Section 4 contains some general remarks about using a principal 
connection to reconstruct an integral curve of a dynamical vector field from one of a reduction 
of it. In Section 5 we describe the two principal connections that can be used in the specific 
reconstruction problem we are concerned with, while in Section 6 we carry out the reduction in 
detail, first in the Abelian case, then in general. In Section 7 we specialize to simple mechanical 
systems, in order to compare our results with those published elsewhere. We conclude the paper 
with a couple of illustrative examples. 

2 Preliminaries 

We will be concerned with Lagrangian systems admitting non-Abelian (that is to say, not nec- 
essarily Abelian) symmetry groups. We begin by explaining what assumptions we make about 
the action of a symmetry group. 

We will suppose that ip M : G x M — > M is a free and proper left action of a connected Lie group 
G on a manifold M. It should be noticed from the outset that this convention differs from the 
one in e.g. j3j [6], but resembles the one taken in e.g. [8| [TO]. 

With such an action, M is a principal fibre bundle with group G; we write M/G for the base 
manifold and tt m : M — > M/G for the projection. We denote by the Lie algebra of G. For any 
£ £ 0) i w ih denote the corresponding fundamental vector field on M, that is, the infinitesimal 
generator of the 1-parameter group V'cxp(^) °f transformations of M. The Lie bracket of two 

fundamental vector fields satisfies [£,fj] = — [£,77] (see e.g. [8]). Since G is connected, a tensor 
field on M is invariant under the action of G if and only if its Lie derivatives by all fundamental 
vector fields vanish. In particular, a vector field X on M is invariant if and only if [£, X] = 
for all £ E 0. We will usually work with a fixed basis for 0, which we denote by {E a }; then for 
X to be invariant it is enough that [E a , X] = 0, a = 1, 2, . . . , dim(0). 
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We suppose that we have at our disposal a principal connection on M. For the most part it will 
be convenient to work with connections in the following way. A connection is a left splitting of 
the short exact sequence 

O^Mxfl^ TM -» (vr M )*T(M/G) -> 

of vector bundles over M; we identify M x pj with the vertical sub-bundle of TM -> M by 
(m, £) i ^ £ | m . Thus we may think of a connection as a type (1, 1) tensor field wonM which is a 
projection map on each tangent space, with image the tangent to the fibre of ir M . The connection 
is principal just when uj is invariant, that is, when C^jj = for all £ € 0. The kernel distribution 
of is the horizontal distribution of the connection. An alternative test for invariance of the 
connection is that its horizontal distribution should be invariant (as a distribution); that is, 
for any horizontal vector field X, [£,X] is also horizontal for all £. We will often refer to a 
connection by the symbol of the corresponding tensor field. 

Let {Xi} be a set of local vector fields on M which are linearly independent, horizontal with 
respect to uj and invariant. Such a set of vector fields consists of the horizontal lifts of a local 
basis of vector fields on M/G, and in particular we may take for the Xi the horizontal lifts of 
coordinate fields on M/G. We then have a local basis {Xi,E a } of vector fields on M. We will 
very often work with such a basis, which we call a standard basis. The Lie brackets of pairs of 
vector fields in a standard basis are 

[Xi,Xj] = RijE a , [Xi,E a ] = 0, [E a ,Eb] = —C^ b E c . 

The Rfj are the components of the curvature of to, regarded as a 0-valued tensor field. The 
second relation simply expresses the invariance of the X; L . In the third expression the C^ b are 
structure constants of Q with respect to the chosen basis. 

It will sometimes be convenient to have also a basis {Xi,E a } that consists entirely of invariant 
vector fields. Let U C M/G be an open set over which M is locally trivial. The projection 
7r M is locally given by projection onto the first factor in U x G — > U , and the (left) action by 
ip™(x,h) = (x,gh). The vector fields on M defined by 

E a : (x,g) ^ ( a d^E a )(x,g) = ^ M {E a (x,e)). 

(where e is the identity of G) are invariant. The relation between the sets {E a } and {E a } can 
be expressed as E a (x,g) = Aa(g)Ef ) (x, g) where (A b a (g)) is the matrix representing ad 9 with 
respect to the basis {E a } of g. In particular, A b a {e) = 5 b a . Since [E a , E b ] = 0, the coefficients A b a 
have the property that E a {Af) = C^Af. 

We revert to consideration of a standard basis. We define the component 1-forms ui a of the 
tensor field ui by ui = uj a E a . Then co a (Xi) = 0, w a (-E^) = 5%. Thus the ui a comprise part of the 
basis of 1-forms dual to the standard basis. We denote by the remaining 1-forms in the dual 
basis. 

Most of the objects of interest, such as the Lagrangian and the corresponding Euler-Lagrange 
field r, live on the tangent manifold of M, which we denote by r : TM — ► M. We recall that 
there are two canonical ways of lifting a vector field, say Z, from M to TM. The first is the 
complete or tangent lift, Z°, whose flow consists of the tangent maps of the flow of Z. The 
second is the vertical lift, Z v , which is tangent to the fibres of r and on the fibre over m coincides 
with the constant vector field Z m . We have Tt{Z c ) = Z while Tr(Z v ) = 0. Moreover, TM is 
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equipped with a canonical type (1, 1) tensor field called the vertical endomorphism and denoted 
by S, which is essentially determined by the facts that S{Z°) = Z N and S{Z V ) = 0. For more 
details on this material, see e.g. [HIS]. The set {Xf , E°, XJ , EY}, consisting of the complete 
and vertical lifts of {Xi,E a }, forms a local basis of vector fields on TM. 

Let {Z a } be a local basis of vector fields on M, and {0 a } the dual basis of 1-forms. These 1-forms 
define fibre-linear functions 6 a on TM, such that for any u E T m M, u = 6 a (u)Z a {m). These 
functions are therefore the components of velocities with respect to the specified vector-field 
basis. We may use these functions as fibre coordinates. Coordinates of this type are sometimes 
called quasi- velocities, and we will use this terminology. In the case of interest we have a standard 
basis {Xi,E a } and its dual {d l ,u) a }; we denote the corresponding quasi- velocities by v % = 
v a = u a . 

We will need to evaluate the actions of the vector fields Xf , E% , XJ and EY on v' 1 and v a . Now 
for any vector field Z and 1-form 9 on M, 

Z C {9) = T^, Z v (8) =t*8(Z). 

Most of the required results are easy to derive from these formulae. The only tricky calculation 
is that of Xf(oj a ), for which we need the Lie derivative of a connection form by a horizontal 
vector field. We have 

(C Xi u a ){E b ) = X i (5Z)-co a ([X i ,E b }) = 0, 
{C Xi u a ){Xj) = -u; a ({X i ,X j ]) = -R? j ; 

in the first we have used the invariance of the horizontal vector fields. In summary, the relevant 
derivatives of the quasi- velocities are 

X<V) = 0, W) = C Xf(v a ) = -R? 3 vi, Xj(v a ) = 0, 
4V) = 0, = 0, 4V) = C&v°, EZ{v b ) = 5 b a . 

Finally, we list some important Lie brackets of the basis vector fields: 

Xf ] = [E a , Xif = 0, [E%,XY] = [E a , X t ] v = 0, 

[ E ai E b] = [Ea,Eb]° = -C c ah E^, [E%,EY} = [E a ,E b ] v = -C c ah E v c . 

3 The generalized Routh equations 

We begin by explaining, in general terms, how we will deal with the Euler-Lagrange equations. 

Consider a manifold M, with local coordinates (x a ), and its tangent bundle r : TM — > M, 
with corresponding local coordinates (x a ,u a ). A Lagrangian L is a function on TM; its Euler- 
Lagrange equations, 

d / dL\ dL () 
di \~dvF J ~ dx® ~ ' 

comprise a system of second-order ordinary differential equations for the extremals; in general 
the second derivatives x a are given implicitly by these equations. We say that L is regular if its 
Hessian with respect to the fibre coordinates, 

d 2 L 

du a du@ ' 
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considered as a symmetric matrix, is everywhere non-singular. When the Lagrangian is regular 
the Euler-Lagrange equations may be solved explicitly for the x a , and so determine a system of 
differential equations of the form x a = f a (x,x). These equations can in turn be thought of as 
defining a vector field T on TM, a second-order differential equation field, namely 

d d 
T = u a ^— + f a - 



dx a du a ' 

we call this the Euler-Lagrange field of L. The Euler-Lagrange equations may be written 

dL \ dL 

0, 



du a J dx c 

and when L is regular these equations, together with the assumption that it is a second-order 
differential equation field, determine T. 

This is essentially how we will deal with the Euler-Lagrange equations throughout: that is, we 
will assume that L is regular and we will work with the Euler-Lagrange field T, and with the 
Euler-Lagrange equations in the form given above. However, we need to be able to express those 
equations in terms of a basis of vector fields on M which is not necessarily of coordinate type. 
It is easy to see that if {Z a } is such a basis then the equations 

T(ZZ{L))-Z*(L) = Q 

are equivalent to the Euler-Lagrange equations. The fact that T is a second-order differential 
equation field means that it takes the form 

r = w a z c a + v a z Y a 

where the w a are the quasi- velocities corresponding to the basis {Z a }. 

We now build in the assumption that L has a symmetry group G, which acts in such a way that 
M is a principal bundle with G as its group, as we described above. We will suppose that the 
Lagrangian is invariant under the induced action of G on TM. This tangent action is defined by 
the collection of transformations ip™ = Tip^ on TM, g £ G. By construction, the fundamental 
vector fields for this induced action are the complete lifts of the fundamental vector fields of the 
action on M; the invariance of the Lagrangian can therefore be characterized by the property 
E%(L) = 0. We have shown in [11] that if L is invariant, then so also is T, which is to say that 

[E^r] = o. 

We choose a principal connection on M, and a basis of vector fields {Xi,E a } adapted to it (a 
standard basis), as described above. Then the Euler-Lagrange equations for L are 

T{XUL))-Xf(L) = 
T{El(L)) - E c a (L) = 0. 

But by assumption E^(L) = 0: it follows immediately that T(E^(L)) = 0. So the functions 
Ea(L) are first integrals, which clearly generalize the momenta conjugate to ignorable coordi- 
nates in the classical Routhian picture. We write p a for E^(L). The Euler-Lagrange field is 
tangent to any submanifold p a = fi a = constant, a = 1,2,... , dim(g), that is, any level set of 
momentum. By a well-known argument (see e.g. |8]), we may regard (x,v) i— > (p a (x,v)) as a 
map from TM to g* , the dual of the Lie algebra Q, and this map is equivariant between the 
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given action of G on TM and the coadjoint action of G on g* (the coadjoint action is defined 
as (£, ad*/z) = (ad g £, fi)). We have 



E%(p b ) = E%E%{L) = [E c a ,El]{L) = -C c ab E-(L) = -C c abPc , 



which expresses this result in our formalism. 

We will also need a less coordinate-dependent version of the Hessian. In fact the Hessian of L 
at w € TM is the symmetric bilinear form g on T m M, m = t(w), given by g(u,v) = u v v v (L), 
where the vertical lifts are to w. We can equally well regard g as a bilinear form on the vertical 
subspace of T W TM, by identifying u and v with their vertical lifts. The components of the 
Hessian g with respect to our standard basis will be denoted as follows: 



We also have g ab = E^(p b ), gi a = XJ (p a ). In general these components are functions on TM, not 
on M, and the Hessian should be regarded as a tensor field along the tangent bundle projection 
r : TM — > M. We will assume throughout that L is regular, which means that g as a whole 
is non-singular. Then V is uniquely determined as a second-order differential equation field on 



We now turn to the consideration of Routh's procedure. We call the function 1Z on TM given 



the Routhian. It generalizes in an obvious way the classical Routhian corresponding to ignorable 
coordinates. The Routhian is invariant: 



whence the result. 

We now consider the Euler-Lagrange equations T(XJ(L)) — Xf(L) = 0. We wish to write these 
equations in terms of the restriction of the Routhian to a level set of momentum, say p a = fi a , 
which we denote by N^. To do so, we need to work in terms of vector fields related to Xf, 
XY and E^ which are tangent to (in general there is no reason to suppose that these vector 
fields themselves have this property, of course). To define the new vector fields we will assume 
that the Lagrangian has an additional regularity property: we will assume that (gab) is non- 
singular. (Note that if the Hessian is everywhere positive-definite then (g a b) is automatically 
non-singular.) Then there are coefficients A 1 -, B\ and C„, uniquely defined, such that 



g(E a , E b ) = g ab , g(X i ,X j ) = g ih g(X h E a ) = g ia = g ai = g(E a , X { ). 



TM. 



by 



K = L- v a p, 



'a 



Eb°(v a pa) = Ci c v c p a - v a C c ba p c = 0, 



(Xf + A\El)( Pa ) 
(Xj + B\El)(p a ) 
(E c a + C b a El)(p c ) 



Xf(p a ) + A b igab = Q 
Xj(p a ) + B b i9ab = 
E c a (p c ) + C b a g bc = 0. 



The vector fields Xf, XJ and E% given by 



X? = X? + A$El 

Aj — Aj + ti i tj a 

E c a = E c a +C b a El 
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are tangent to each level set N^. (The notation is not meant to imply that the barred vector 
fields are actually complete or vertical lifts.) We will need to know the coefficients explicitly 
only in the case of Bf and C h a : in fact 

Bf = -9 ah 9ib and & a = g bc C d acPd . 

This is all carried out under the assumption that (g a b) is non-singular. One has to make such 
an assumption in the classical case in order to be able to solve the equations dL/d9 a = \i a for 
the 9. In the general case the non-singularity of (g a b) is the condition for the level set to 
be regular, i.e. to define a submanifold of TM of codimension dim(g). The vector fields E% are 
transverse to all regular level sets, and the barred vector fields span the level sets. Thus on any 
regular level set the bracket of any two of the barred vector fields is a linear combination of vector 
fields of the same form. We want in particular to observe that this implies that [E^,XJ\ = 0. 
It is not difficult to see, using the known facts about the brackets of the unbarred vector fields, 
that this bracket is of the form P a E^; this must satisfy P a E^(pb) = P a g a b = 0, whence by 
the regularity assumption P a = 0. In fact by similar arguments the brackets of the barred 
vector fields just reproduce those of their unbarred counterparts, except that [X^,XJ] = 0. In 
particular, [E°,E£] = —C^ b E^. The E% therefore form an anti-representation of g, acting on 
the level set (just as the E a do on M). 

We return to the expression of the Euler-Lagrange equations in terms of the Routhian. We will 
need to evaluate the actions of Xf and X? on v a . Using the formulae in Section 2 we find that 

X?(v a ) = {Xf + A\El){v a ) = -Rt^+A1 
X?(v a ) = {X- + BlEl){v a ) = B?. 

We now set things up so that we can restrict to the submanifold easily. We have 

X?{L) = Xf{L)-A*El{L) 

= Xf{L - v a p a ) + (-i?»V + A?)p a + v a Xf (p ) - Atp a 

= Xf{K)-p a R a ijV i- 

X?(L) = XUL)-BtK{L) 

= X^ (L - v a Pa ) + Btpa + v a X( (p ) - Btpa 

= xnn). 

But T(XJ(L)) - Xf{L) = 0, and T is tangent to the submanifold iV M ; thus if we denote by W 1 
the restriction of the Routhian to the submanifold (where it becomes L — t> a /i a ) we have 

T(XUnn) - X?{K») = -ii a R%vK 

On the other hand, if T is a second-order differential equation field such that T(E^(L)) = 
and the above equation holds for all fj, a then T satisfies the Euler-Lagrange equations for the 
invariant Lagrangian L. 

We will refer to these equations as the generalized Routh equations. 

Neither TV 1 nor T is ^-invariant. They will however be invariant under those vector fields £ c , 
£ G g, which happen to be tangent to the level set iV^. These are the vector fields for which 
£, a E2 = £ a E2, or ^ a C^ b fi c = 0. We will return to this issue in later sections. 
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Note that since T satisfies T(p a ) = it may be expressed in the form 

v = v i xf + r i i I v + /E a c . 

If the matrix- valued function XJX^ijV) is non-singular, the reduced Euler-Lagrange equations 
above will determine the coefficients P. We show now that this is always the case, under the 
assumptions made earlier. 

Recall that X? = X? + BfEJ is determined by the condition that XJ(p a ) = 0, and that 
therefore Bf = —g ab gi b . We may regard Xj + BfE a as a vector field along the tangent bundle 
projection, and XJ really is the vertical lift of this vector field; we will accordingly denote it by 
Xi. Then 

g(Xi,E a ) = g{Xi, E a ) + B b g{E b , E a ) = g ia + B\g ah = 0; 

thus the Xi span the orthogonal complement to the space spanned by the E a with respect to the 
Hessian of L. That is to say, the tangent space to a regular level set of momentum at any point 
u € TM intersects the tangent space to the fibre of TM — > M at u in the subspace orthogonal 
with respect to g u to the span of the EJ. Moreover, 

g(Xi,Xj) = g l] +Btg aj + B^g ia + BtB b J g ab 
= gij - 2g ab g ia g jb + g ac g hd g ic gjdgab 

= 9ij — g ab gia,gjb- 

So this is the expression for the restriction of the Hessian of L to the subspace orthogonal to 
that spanned by the E a . 

Now recall that XJ(1Z) = X?(L). Thus 

XJXJ(TZ) = (XY - g ah g lb El)X](L) = 9ij - g ab g ib g aj = g(X u Xj)- 

That is, the 'Hessian' of 1Z (i.e. XjX^ijV)) is just the restriction of the Hessian of L to the 
subspace orthogonal to that spanned by the E a . It follows that the bilinear form with compo- 
nents g^ = XjXjilV) must be non-singular. For suppose that there is some vector w J such 
that (jijwi = 0; then g(Xi,w^ Xj) = by assumption, and g(E a ,w^ Xj) = by orthogonality - 
but then w J Xj = since g is assumed to be non-singular. 

The sense in which the generalized Routh equations are 'reduced' Euler-Lagrange equations is 
that (in principle at least) we can reduce the number of variables by using the equations p a = \i a 
to eliminate the quasi-velocities v a . However, these variables appear explicitly in the expression 
for r, so it may be considered desirable to rearrange the generalized Routh equations so that 
they no longer appear. This can be done by changing the basis of vector fields on the level set 
of momentum, as follows. The change is suggested by the fact that, notation notwithstanding, 
S(X?) ^ XJ (where S is the vertical endomorphism). Let us, however, set 

V"C v"C i rja i^C . 

A i — X i + B i &a ■ 

then since S vanishes on vertical lifts, 

S(X?) = S(X° + B?E c a ) = X? + B1K = XJ. 

We write 

r = jxf + rxv, 
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so that 

r = r + (v i Bf + v a )EC. 

We will examine the contribution of the term involving E% in T to the generalized Routh 
equations. First we determine E%(Jl). Since E^{p b ) = 0, 

E c a (1l) = E c a (L - v b Pb ) = C b p b - E c a (v b )p b = C b aPb - C b acPb v c - C c J b p b = -C b aJ?b v c . 

It follows that 

E c a (x-(ii)) = xumn)) = -XUC b ac p b v c ) = -C b aJ)b Bf. 
So setting T = T + (v*Bf + v a )E% we have 

T{XJ{^)) ~ Xf{K») = r (X i v (^)) + (v j B% + v a )Ef(X?(R,»)) - Xf{K») + BfEf(TZ^) 

= T (X?(Kn) - Xf{K») - {v*Bf + v a )C h ^ b Bt - B?C b ac p b v c 
= r (X 4 v (^)) - Xf{W) + viB*BiC b ac ^ 

and the generalized Routh equations become 

r (X 4 v (7^)) - xf{Rr) = -p a {R% + B b B^ c y. 

We may say that among the vector fields tangent to a level set of momentum it is Xf + BfE^ , 
not Xf, that really plays the role of the complete lift of Xi. Be aware, however, that unless the 
symmetry group is Abelian, To cannot be identified with a vector field on T{M/G). We will 
return to this matter at the end of Section [5j 

To end this section we give a coordinate expression for the generalized Routh equations in their 
original form. For this purpose we take coordinates (x l ) on M/G, and coordinates (x l ,9 a ) on 
M such that the 9 a are fibre coordinates; then (x l , 6 a ,v' 1 ) are coordinates on iV„, which is to say 
that can be locally identified with M Xm/g T(M/G). We may write 

Y - A _ \*JL P - — 

'~ dx* '89-' a ~ a d9 b 

for suitable functions A" and K b on M. (We should note that the K b are components of a 
non-singular matrix at each point; moreover, the invariance property of the Xi can be expressed 
in terms of the coefficients A" and K b ; but we will not actually need either of these facts here.) 
From the formulae for the action of complete and vertical lifts on quasi-velocities given at the 
end of Section 2 we see that 

Xf («*) = Eftf) = 0, X?{vi) = 6> i . 
Thus in terms of x l , 6 a and v 1 we can write 

vc = _^ A a -^- X v = — F c = K h - 

* dx* i de a ' 1 dv v a a de b ' 

It is necessary to be a little careful: the coordinate vector field expressions are ambiguous, 
since they can refer either to coordinates on TM or on N^. We emphasise that it is the latter 
interpretation that is intended here. In view of the possibilities of confusion it will be useful to 
have an explicit notation for the injection N„ — > TM: we denote it by l. The non-singularity 
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of (gab) ensures that, at least locally, we can rewrite the relation p a = fj, a for the injection 
i : — ► TM in the form v a = i a (x\ 8 a ,v % ), for certain functions i a of the specified variables. 

The restriction of the Euler-Lagrange field T to iV M is 

r = i a M% + jxf + (r* o L )xy 

= ^w^'(l>-^w) + ^^ a 



dv i 

d9 a ' " dx { ' ^ ' "'dv 1. 



the equations for its integral curves are 



x = v , 

9 a = L b (x,9,v)K%(x,6)-v i h a i (x,9). 
These can be considered as a coupled set of first- and second-order equations, 

J? = P(x,0,x), 

\ e a = i b (x,e,x)K^(x,e) -x i Af(x,e). 

With regard to the second of these equations, we point out that the expression for the velocity 
variables in terms of the quasi- velocities v l and v a is just 9 a = v b K£ — v i Af. What turns these 
identities into genuine differential equations is, in particular, substitution for the v a in terms of 
the other variables via the functions i a — or in other words, restriction to N^. 

The functions P may be determined from the generalized Routh equations. These may be 
expressed as 

dt \ dv* J dx* Ma ij * de a ' 

In the light of the earlier remarks about the interpretation of coordinate vector fields, we point 
out that substitution for v a in terms of the other variables in this equation must be carried out 
before the partial derivatives are calculated. 



4 The reconstruction method 

We have seen in Section 3 that Routh's technique consists in restricting the Euler-Lagrange 
equations to a level set of momentum N^. This procedure takes partial, but not necessarily 
complete, account of the action of the symmetry group G. To make further progress we must 
examine the residual action of G on N^. 

As we mentioned before, the momentum map is equivariant between the induced action of G on 
TM and the coadjoint action of G on Q* . The submanifold is therefore invariant under the 
isotropy group G M = {g € G | ad* fj, = fj,} of fx. The algebra of G M consists of those £ G 
such that t^C^bHc = 0; this is the necessary and sufficient condition for £ c to be tangent to N^. 

Note that any geometric object we know to be G-invariant is automatically G^-invariant. 

The manifold is a principal fibre bundle with group G^; we will denote its base by N^/G^. 
The restriction of the Euler-Lagrange field V to is G M -invariant, and as a consequence it 
projects onto a vector field T on N^/G^. 
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The task now is to examine the relationship between f and T. There are two aspects: the 
formulation of the differential equations represented by T; and the reconstruction of integral 
curves of T from integral curves of T (supposing that we have solved those equations). 

Our methods of attack on these problems will be based on those we developed in our papers [3llllj 
and are similar to (but different from) the ones that were adopted in e.g. [7]. These in turn were 
based on the following well-known method for reconstructing integral curves of an invariant 
vector field from reduced data. Let ir : N — > B be a principal fibre bundle with group G. 
Any invariant vector field T on N defines a -zr-related reduced vector field f on B: due to the 
invariance of T, the relation Tir(T(n)) = r(7r(ra)) is independent of the choice of n E N within 
the equivalence class of 7r(n) G B. Given a principal connection O, an integral curve v(t) of T 
can be reconstructed from an integral curve v(t) of T as follows. Let v H (t) be a horizontal lift 
of v(t) with respect to Q (that is, a curve in ./V over v such that £l(v H ) = 0) and let g(t) be the 
solution in G of the equation 

m*)) = ^owt))) 

where $ is the Maurer-Cartan form of G. (We use here the fact that given any curve in 
Q, the Lie algebra of G, there is a unique curve g(t) in G which satisfies $(g(i)) = £(i) and 
g(0) = e; g(t) is sometimes called the development of £(t) into G, see for example [H].) Then 
v(t) = tpg^v H (t) is an integral curve of V. 

In the following sections we define two principal connections on Np, we determine T and we 
identify for both connections the vertical part of T, necessary for the reconstruction method 
above. 

5 Two principal connections on a level set of momentum 

A principal connection on — > N^/G^ is by definition a left splitting of the short exact 
sequence 

-» Np x 0/i - TNp - AV x^ /Gm T(NjGp) - 0; 

all spaces in the above sequence should be interpreted as bundles over N^. We think of Q as a 
type (1, 1) tensor field on which is pointwise a projection operator with image the tangent 
space to the fibre, and which is invariant under G^. 

The first connection we define uses the Hessian of L to determine its horizontal distribution, 
and is therefore analogous to the mechanical connection of a simple system; we denote it by O m . 

Recall that we interpret the Hessian g of L as a tensor field along r. In particular, its components 
with respect to the standard basis {Xi,E a } are functions on TM. We will say that a vector 
field W on is horizontal for f2 m if 

g(i,T*W) = 0, V£e 0M , 

where t*W is the projection of a vector field W on TM to a vector field along r : TM — > M. 
The definition makes sense only if we assume that the restriction of g to x is non-singular, 
as we do from now on. 

In [11] we have shown that if the Lagrangian is invariant then so is g, in the sense that 

Cxg = 0, V£ € 0. 
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Here, for a vector field Z on M, £z stands for an operator acting on tensor fields along r that 
has all the properties of a Lie derivative operator, and in particular, when applied to a function 
/ on TM and a vector field X along r gives 

c f-7°(f\ c x- (yP dxa + dz \n dxa x? 92 ^ 8 

£zt-Z{f), £ Z X - \Z^— T + — — u'— T -X^- 



dx? dx^i dul 3 dx? / dx a 



where 

Z = Z°{x)-£-, X = X°(x,u) ° 



dx a ' v ' ' dx a ' 

Note that if X is a basic vector field along r (i.e. a vector field on M), then CzX = [Z,X]. 
Furthermore, for any vector field W on TM we have 

£ z (t*W) = t*[Z c ,W}. 



To show that the connection is principal we need only to show that if W is horizontal so also is 
[|°, W] for all £ G 0„. But for all £, V € 0„, 

5(n[| c , W],»7) = 5 (£ f (r^),?7) = -5(r,W,^) = -g(r*W, [Ifj]) = g(nW, [^?]) = 0, 

using the properties of the generalized Lie derivative and the invariance of g. 

As was mentioned before, {Xf, E%, X^} is a basis of vector fields on N^. Suppose now that the 
basis {E a } = {Ea,E } of g is chosen so that {Ea} is a basis of g^. Then C Ab fi c = 0, and on 
Nn we get for the corresponding fundamental vector fields 

E A = E A + 9 hc C\ c Vd&b = E% 

All E^ are therefore tangent to N„, as required. These vector fields span exactly the vertical 
space of — > N^/G^ which we have identified with x g^. Vector fields of this form are 
infinitesimal generators of the G M -action on N^. 

If (G AB ) is the inverse of the matrix (gAs) (and not the (A, S)-component of (g ab )), then the 
vector fields 

pH _ BC ^AB pC _ 771C x 6 F c 

v-h -^-c r<AB „ So y-c po 

Aj — Aj - G gAi&B — A i - l i 

together with X? , are horizontal. (As was the case with the notations E% etc., the notation for 
the horizontal fields is not meant to imply that E™ etc. are actually horizontal lifts.) The action 
of S7 m is simply 

n™(E c A ) = e% ^{ed = o, n™(xf ) = o, n m (x?) = o, 

and since the arguments form a basis of vector fields on these equations specify £l m explicitly. 
We will call f2 m the mechanical connection on iV„. 

The vector fields Xf = Xf + BfEf introduced earlier are also horizontal; they can be expressed 
as Xf = Xf - g ab g H El. 

The vector fields Xf are not horizontal with respect to £l m . However, it is possible to identify 
a second principal connection Q N f on for which these vector fields are horizontal. We will 
identify in two steps. 
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It seems natural to split the basis {X^,E^,X^} into a 'vertical' part {E%} and a 'horizontal' 
part {X?,Xy}. To see that it does indeed make sense to do so it is sufficient to observe that the 
distributions spanned by {E%} and {Xf ,X^}, respectively, are unchanged when the bases {E a } 
of Q and {Xi} of ^-horizontal vector fields on M are replaced by different ones. Under a change 
of basis for the {E% } are simply replaced by constant linear combinations of themselves, so 
their span is clearly unchanged. On the other hand, if we set Yi = A\Xj then = AjXJ 

and YJ° = A\X^ + A^XJ (where A\ is the total derivative of A?, not that it matters), so the 
distributions spanned by {X^,Xf} and {Y^ C ,Y^ V } are the same. 

So we can indeed characterize a connection in this way, but it is not a connection on Nu — > 
Nn/Gu- In fact this construction defines a connection on the bundle with projection Nn — > 
T{M/G) (the restriction of Tir M : TM -> T(M/G) to N^), i.e. a splitting of the short exact 
sequence 

- N„ x g -> TN U -» JV„ x T(M/G) T(T(M/G)) -> 0. 

(Recall that the vector fields E°, which span the vertical space of the projection Nu — > T(M/G), 
form an anti-representation of acting on the level set iV M .) The construction just described is a 
version of the so-called vertical lift of a connection on a principal bundle (here u) to its tangent 
bundle (this is described more fully in [3]); accordingly we denote the corresponding type (1,1) 
tensor field by f2 v , and we have 

n v (E°) = E c a , tt v (xf ) = o, n v (x?) = o. 

Evidently (ft v ) 2 = n v . We show now that £gcf} v = for all A. Firstly, note that 

[pc pci rpc ttici r^> i? c 

L-^Aj-^oJ - \E A ,tj a \ - -L Aa ti b , 

so that 

(£^ v )(E C a) = [E%^m\ - n v [E c A ,E c a ] = [E%E°] + C b Aa E£ = 0. 
Moreover, since [E%X°] = 0, 

(£^cO v )(Xf ) = [iS,o v (xf)] - O v [^,Xf] = 0, 

and similarly for XJ . 

The relation between £l v and w may be described more easily if we momentarily break our 
convention by specifying connections by their forms rather than by the tensors corresponding 
to their splittings: it is easily checked that 

n v (Z v ) = u;(T(TOi)Z v ), Z v € TN^. 

This equation has to be read as one between elements of 0, obtained by identifying the vertical 
subbundle of TNu with N^ x g and the vertical subbundle of TN with M x g, or if you will by 
projection onto 0. 

The vertical space N^ x g^ of the connection Q, N i* we are looking for is only a subbundle of the 
vertical space N^ x g of the connection Q v . So in a second step we need to identify a connection 
for the following sequence of trivial vector bundles: 

-> iV M x 0M -» iV M x -► N^ x 0/ 0M -> 0. 
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For this connection we can simply take the restriction of the mechanical connection f2 m defined 
earlier to the submanifold x g. The connection Q N i* is then simply f2 m ofi v (see the diagram 
below) . 

N» x M - N^xg ^ x 



N„x T(M/G) T(T(M/G)) - N„x T(M/G) T(T(M/G)) 

By construction ft v of} m = Q m , so fi^f = fi m of] v satisfies (O^) 2 = as it should. We have 

n N " (e c a ) = e% n N » {ed = o, si N " (x?) = o, n N » (A7) = o. 

The tensor field Q, N ^ is G^-invariant since both of the tensors of which it is composed are 
G^-invariant; Q N » therefore defines a principal G^-connection. 

Note that to define the mechanical connection we do not need a principal connection to on 
M — > M/G (though we may use one in calculations). If such a connection is available then we 
can use either O m or Q N i* for the reconstruction method. 

The connection Q, N » is clearly different from Q m in general. We can also decompose f2 m into 
two connections, in accordance with the short exact sequences in the diagram. The splitting 
Qq of the middle vertical line, similar to the connection Q Y of £l Nfl , can be defined by saying 
that a vector field W is horizontal if g(!;,T*W) = for all £ £ g (not just for £ £ g ). For this 
connection we have 

n v (E c a ) = E c a , nytf?) = o, n v (x?) = o, 

where the vector fields Xf are exactly those that we have encountered in Section 2. 

To end this section we consider the decomposition of the restriction of the Euler-Lagrange field 
r to into its vertical and horizontal parts with respect to the two connections. 

Let us introduce coordinates (x l , 9 a ) on M such that the orbits of G, or in other words the fibres 
of M — > M/G, are given by x % = constant; the x % may therefore be regarded as coordinates on 
M/G. As before, we will use as fibre coordinates the quasi-velocities (v l ,v a ) with respect to 
the standard basis {Xi,E a }. The non-singularity of (g a b) ensures that, at least locally, we can 
rewrite the relation p a = fj, a for the injection i : — > TM in the form v a = L a (x l ,9 a ,v l ), for 
certain functions i a of the specified variables. The restriction of the Euler-Lagrange field to 
is 

r = i a E c a + v'x? + (r* o L )xy 

= {L A + rinE c A + L a Ei + v*x? + {ro i )xj 

= (t A + T A i a + Tfv^Ei + i a El + jxf + (P o l)XJ. 
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The vertical part of T with respect to the mechanical connection fi m is (l a + T A i a + Tfv 1 )^^, 
and with respect to the vertical lift connection Q N f it is (l a + T A L a )E^. 

Note that neither of the current decompositions of T coincides with the one we had towards the 
end of Section 3, which we should now write T = (i a + Bfv l )E^ + IV The reason is that this 
last decomposition is only partial, in the sense that the vector field Tq is the horizontal part of 
r with respect to the connection Qq ; it is the horizontal lift of a section of the pullback bundle 
Nfj, Xt(m/G) T(T(M/G)), not a vector field on T(M/G), and this section is only a part of the 
data required for the reconstruction method. 

6 The reduced vector field 

A principal connection is all we need to reconstruct integral curves of an invariant vector field 
from those of its reduced vector field. We next examine the latter. 

6.1 The Abelian case 

Before embarking on the more general case, it is instructive to see what happens if the symmetry 
group G happens to be Abelian, i.e. when C^ b = 0. Then as we pointed out earlier for the case 
of a simple mechanical system with Abelian symmetry group, 0^ = and any level set p a = (i a 
is invariant under the whole group G. In fact, under the assumption that p a = \i a can be solved 
locally in the form v a = t a , N^/G can be interpreted as T(M/G), with coordinates (x\t> 1 ), 
where the x l are coordinates on M/G and the v l the corresponding fibre coordinates (no longer 
quasi- velocities) . In this case there are no l E a '-vectors and E% = E% for all a. 

The restriction of the Euler-Lagrange field to N^, given here by 



is now also G-invariant. As a consequence, the coefficients P o i do not depend on the group 
coordinates 6 a but only on the coordinates of T(M/G). In fact the vector fields v l Xf + 

(P o t)xy (the n N " -horizontal part of T) and v*X? + (P o l)XJ (the ^'"-horizontal part of 
r) both reduce to the same vector field on N^/G, which in this case is exactly a second-order 
differential equation field on T(M/G), 



The integral curves of this reduced vector field are the solutions of the equations x l = P(x,x) 
(with v % = x 1 ) and, from the introduction, we know that these are equivalent to the equations 



6.2 The non- Abelian case 

In the general case of a non-Abelian symmetry group we should not expect that the equations 
for x l will be completely decoupled from all coordinates 6 a . Indeed, in that case the vector field 



r = i a E° + jxf + (p o t )xy 

= {t a -g ab guv i )E c a +v i Xf + {Y i o L 
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r reduces to a vector field on N^/G^. This manifold can locally be identified with M/G M x 
T{M/G), so the equations for the integral curves of the reduced vector field will depend also on 
the coordinates of M/G^. 

To give a local expression of the reduced vector field T we need to introduce a basis for 
X {N n / G fj) . The bracket relations [E%Xf] = and [E%XJ] = show that Xf and Xf 
are G^-invariant vector fields on N^; they project therefore onto vector fields Xf and Xf 
on N^/G^. The invariance of the Hessian g amounts for its coefficients to E°(gij) = 0, 
K (9bc) + C d ab g cd + C d ac g M = 0, and E° (g ib ) + C c ah g tc = 0. From this 

E c A (?f) = c% c rf, 

f?C / r v*B\ /~iB <Y"C r-ift y^B s~iB 

^A\ L a) — L/ AC l a- L/ Aa l P ~ L/ Aa 

(where we have taken into account the fact that in the current basis C AB = 0). It is now easy 
to see that the vector fields Xf = Xf — Y A E A are also G^-invariant. In fact, since they differ 
from Xf only in a part that is vertical with respect to the bundle projection — > N^/G^, 
they project onto the same vector fields Xf on N^/G^. 

The vector fields Ef are not invariant: in fact [E A ,E^] = —C Aa E^. To obtain a complete basis 
for X(Nn/Ga), we need to replace the vector fields {Ef,} by G^-invariant vector fields. To do so, 
we will consider the G- invariant vector fields E a = A^Ef, on M that we introduced in Section 2. 
Let (Ap) be the coefficients we find in the relation E a = A^Ep + A^Eb- The vector fields 

k = Am 

are tangent to the level set 7V„ and horizontal. Given that C AB = 0, it easily follows from the 
relation Ea(A^) = C Al Al that these vector fields are G^-invariant: 

[E%Ef] = E A {Ai)E^ - AlC\pE* = 0. 

They project therefore onto vector fields Ef on N^/G^. To conclude, the set {Xf,Xf,E^} 
defines the basis for X(N fl /G fl ) we were looking for. 

We denote by (A^) the matrix inverse to (Af) and set $ A = i A + T A i a , <5> A = l a + T A i a + T^V 
and ^ a = A a p ft '. Then T takes the form 

r = ^E^ + ^El + tfXf + (P oi)Xf, 

where, as before, the first term is the vertical part of V with respect to the vertical lift connection 
Q, N n (i n the first place) and the mechanical connection Q m (in the second). Obviously v i and 
(P o l) are G^-invariant functions. To see that ^> a is also G M -invariant, recall that C AB = 
and E c A {tP) = C^C, and observe that E%(A$)A} = -A P S E C A {A}) = -A^A%C\ a = -CT M . 
Therefore 

E C A (V*) = E c A (A$y + Affitf) = -A a 5 C 5 A / + A$C^ = 0. 

We conclude that t>% (P o i) and \E ,a can all be regarded as functions on N^/G^. The horizontal 
part of T, for both connections, can thus be interpreted as the horizontal lift of the reduced 
vector field 

f = ^ a Ef + v l Xf + (P o L )Xf 



17 



on N^/G^. 

For completeness, we point out that it follows from the relations E^(T B ) = C B C T^ — C^ a T B - 



Cf a and E%l b ) = C% c i c + C^P that the coefficients $ A and $ A satisfy 

E C A (<S> B ) = C% c t> c , E%($g) = C% c <t>g. 

This shows that they can be interpreted as the coefficients of g^-valued functions $ and <J> m 
on satisfying $ o ip g N>1 = &d g <3? for g £ (and similarly for <l? m ), where ip Nli denotes the 
G^-action on (see [3]). 

We now give a coordinate expression for the reduced vector field. From here on we will use 
coordinates (9 a ) = (9 A , 9 a ) such that the fibres of G -> G/G^ are given by 9 a = constant. With 
this assumption, there are functions on M such that 



d9 B> d6 B d9@ 



We also introduce the functions A b for which 



X = — - \ b — 

1 dx i i d9 b ' 



as before. 



By interpreting N^/G^ locally as MjG^ x M / G T(M/G), we see that a point of N^/G^ has coor- 
dinates (x l , 9 a ,v l ). Because of their G^-invariance, the functions Pot and *$> a are independent 
of the variables 9 A . Let tt n ^ be the projection Nn —> N^/G^; then for any invariant function 
F on Nfj, there is a function / on N^/G^ such that F = / o Then for all invariant vector 
fields X on iV„, and their reductions X to vector fields on N^/G^, we have 



= X(f) O 7T 

We will apply this property to the vector fields E^, Xf and Xf and the invariant functions x % 
v % and 9 a . Keeping in mind that for any vector field Z, function / and 1-form 9 on M, 



Z c (r*f) = r*Z(f), Z v (r*f) = 0, Z c (9) = £ z 9, Z v (9) = t*9(Z). 

where 9 stands for the fibre- linear function on TM defined by the 1-form 9, and r is the tangent 
projection TM — > M, we find that 

X?( X j) = 6i Xf(^) = -Af, Xf(^) = 0, 

El{xl) = 0, E«(9?) = AZK?, = 0, 

from which it follows immediately that 

and 

f = (J*k% - u*Af + + (r j o t )^. 
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The equations that determine the integral curves i>(t) = (x l (t), a (t), v l {t)) of the reduced vector 
field f 1 are therefore the coupled set 

X* = P O L, 

One can easily convince oneself that the right hand side of the equation for 9 a is indeed inde- 
pendent of the variables 9 A : by considering the coefficients of d/dd a in [E A ,Ep\ = —C A ^E C we 

find that E A (Kf) = -CP Ap K^ since also E^{iP) = C%p and E A (A?) = 0, it follows easily 

that E%(iPK% - v { Af) = 0, as claimed. 

The functions P o i on the right-hand side of the equation for the 'x % can be determined from 
the generalized Routh equations of Section 3, 

T{XUW))-Xf{K») = -ii a Rt^. 

Since W 1 = TZ o t is G^-invariant, so also are X?(K^) and Xf{TZ^). Recall that the Rfy are 
functions on M, determined by = RfjE a . Thus since [E a ,Xi] = 0, 

(E a (Rij) ~ RijCac) E b = 0. 

But C Ac Hb = 0, and so 

E c A {R\^ b ) = E A {R%)ii b = Rtf\ c[ i h = 0. 

It follows that the term fiaRfjV 3 is G^-invariant. The generalized Routh equations therefore 
pass to the quotient N^/G^, and take the reduced form 

r(xnnn)-x?(nn = -»aR?y. 

Following [10], we will call these reduced equations the Lagrange-Routh equations. Under the 
regularity assumptions we have adopted throughout, the function- valued matrix (X^ XJ (TZ^)) is 
non-singular and the coefficients Pot, now interpreted as functions on N^/G^, can be determined 
from the Lagrange-Routh equations. In the current coordinate system the equations become 

dt \ cv J dx* ~ Ma » 1 de<* ' 

Given a reduced solution v(t) = (x l (t) , 0° (t) , v l (t)) £ N^/G^, we can apply the method of 
reconstruction using either one of the connections O m and D, * to recover a complete solution 
v(t) = (x i (t), 9 A (t), 6 a (t), v l (t)) € of the Lagrangian system. The examples discussed in 
Section [8] will make it clear how this method works in practice. 

7 Simple mechanical systems 

In this section we reconcile our results with those for the case of a simple mechanical system to 
be found elsewhere in the literature. 

A simple mechanical system is one whose Lagrangian is of the form L = T — V where T is 
a kinetic energy function, defined by a Riemannian metric g on M, and V is a function on 
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M, the potential energy. The symmetry group G consists of those isometries of g which leave 
V invariant. We define a connection on M — > M/G by taking for horizontal subspaces the 
orthogonal complements to the tangent spaces to the fibres; it is this connection that is usually 
called the mechanical connection (for a simple mechanical system). We write g a b = g(E a ,Eb), 
gij = g(Xi,Xj); by assumption, g(E a ,Xi) = 0. Then in terms of quasi-velocities, 

L(m,v) = \ ( y g i j{m)v l v j + g ah (m)v a v h ^j - V(rn), 

and E^E^(L) = g a b etc., so the notation is consistent with what has gone before. Note that since 
we assume that g is Riemannian and therefore positive-definite it is automatic that L is regular 
and that the matrices (g a b( m )) an d (gij(m)) are both non-singular for all m; in particular we 
don't need to make the separate assumption that (gab) is non-singular. Considered as defining 
a map M — > g* Q* , (g a b) is called the locked inertia tensor. The isometry condition gives 

E a (9bc) + C d ah g cd + C d a C g hd = 0, E a ( gij ) = 0. 

The first of these is the differential version of the equivariance property of the locked inertia 
tensor with respect to the action of G on M and the coadjoint action of G on g* g*. The 
second shows that gij may be considered as a function on M/G. 

The momentum is given simply by p a (m,v) = g a b( m )v b - On any level set N^, where p a = fi a , 
we can solve explicitly for the v a to obtain v a = g ab ^b- 

The Routhian is given by 

11 = L- Pa v a = Igijjv* - \g a bv a v h - V, 
and on restriction to Nn we obtain 

K" = \g i3 v L vi - (V + ±<? a W fe ) • 

The quantity V + \g ab ^ a ^b is the so-called amended potential [15] and the term = \g ah \ia^b 
is called the 'amendment' in [TO]. Both functions on M are G^-invariant: one easily verifies that 
E a (C^) = g^Cab/^c^d, so in particular for a = A we get Ea(C^) = 0. 

Note that by the choice of connection Bf = 0; we have Xf = Xf , and the generalized Routh 
equations are 

T (X?(K>*)) - Xf{TW) = 

This equation is the analog in our framework of the one in Corollary III. 8 of [TO]. We have 
shown in the previous section that it reduces to the Lagrange-Routh equations 

f(X?(W))-X?(W) = -» a RZ j vi, 

which for consistency should be supplemented by the equation that determines the variables 
6 a . As we pointed out earlier, the latter is actually just the expression for genuine velocity 
components 9 a in terms of quasi-velocities, supplemented by the constraint v a = L a which for a 
simple mechanical system takes the form i a = g aa ji a - 

We can split the reduced Routhian W 1 in a Lagrangian part £ = \gijV % v^ — V and the reduced 
amendment (2^. Since the quasi-velocities v a do not appear in the expression of £, it can formally 
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be interpreted as a function on T(Q/G). The reduced amendment is a function on Q/G^. We 
can now rewrite the Lagrange-Routh equations as 

f (*y(£)) - Xf (£) = -/^V + 

in coordinates 

s(£)-§--^ + 

This equation is only one out of two equations that appear in Theorem III. 14 in [10], the theo- 
rem that states the reduced equations obtained by following a variational approach to Routh's 
procedure. We leave it to the reader to verify that the second equation, in its form (III. 37), is 
in fact 

v a C^a = g ab v b C a ^ a . 

Since v a = g ab ^b-, this is obviously an identity from the current point of view; it certainly 
cannot be used to determine 9 a in terms of the other variables, and without this information 
the equations are incomplete. In this respect, therefore, our reduction results are an improvement 
on those in 10 1. 



Let us now check, in the case where the configuration space M is of the form S x G, for an 
Abelian symmetry group (C^ b = 0), and a Lagrangian of the form 

L{x,9,x,9) = \k ij {x)x i x j + k ia (x)± i 9 a + \k ah {x)9 a 9 b -V(x), 

that the reduced equations above coincide with those in the introduction. We set 

W - K b — 

a ~ a d9 b 

where the are independent of the 9 a since we are dealing with the Abelian case. In general, 
horizontal vector fields take the form 

X = — -A a — 

1 dx i 1 89 a ' 

The quasi-velocities adapted to the connection are therefore given, as before, by v l = x % and 
K a y b = 0a + A.fx\ 

Given that in this case E a (gb c ) = and E a {gij) = 0, all coefficients of the metric can be 
interpreted as functions on M/G = S, and they depend only on the variables x 1 . The use of the 
mechanical connection entails that g a { = 0. When expressed in terms of the coordinates (x l ,9 ), 
this property fixes the connection coefficients to be of the form A" = k ab kib and the remaining 
coefficients of the metric to be g a b = k ca \K c a K^ and = kij — k ab ki a kji). The expression for W" 
given in the introduction now easily follows. Since TV 1 is a function only of x % and v l = x l we 
get 

X J(K») = ^r- and Xf (TV) = 

1 v ; dx l 1 y ! dx l 

Moreover, 

d 

[Xi,Xj] = BijQQa- = RijEa, 

and likewise /j, a = E^(L) = K b TTb- So ^aRfj = ^aBfj and the equation from the introduction 
follows. 
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We now return to the general case (of a simple mechanical system) and consider the reconstruc- 
tion process. 

We continue to use the mechanical connection on M. Since, in the basis that is adapted to this 
connection, gi a = 0, and therefore Xf = Xf , the two connections Q m and Q Nf * coincide. We 
denote the common connection on by fi. 

Let v(t) be a curve in N^/G^ which is an integral curve of T, and v H a horizontal lift of v to 
(horizontal with respect to Q). The reconstruction equation is 

Hjif = n(r o t) H ) 

(where i? here is the Maurer-Cartan form of G^); this is (at each point on the curve v) an 
equation between vertical vectors on N^, but can and should be thought of as an equation on 
g^. It determines a curve g{t) in G^ such that 

t » ^v«(t) 

is an integral curve of T in Nu) again, ip Nli is the action of G^ on N^. So far, this works for an 
arbitrary Lagrangian. 

Now as we showed earlier in general, the vertical part of T with respect to the vertical lift 
connection is (l a + T A L a )E% Thus in the case at hand 

O(r) = {i A + TiOE c A 

= (i A + G AB g Ba L a )E c A 

= (g A y a + G AB g Ba g a y a )E% 

= (g A >a + G AB (5 B -g BC g Ca )Vam 
= (/> a + G A V-2 A >a)^ 
= G AB hbE% 

The first point to note is that the coefficient G AB [i B appearing on the right-hand side of the 
final equation above is a function on M, so that in the right-hand side of the reconstruction 
equation the argument v H can be replaced by its projection into M, which is roto v H . 

Next, we interpret G AB \x B in terms of the locked inertia tensor. Recall that the locked inertia 
tensor atm£M has components gab(™). As is the usual practice we consider the locked inertia 
tensor as a non-singular symmetric linear map I{m) : g — > g*. Now let j be the injection 
g M — > g: then [i B are the components of j*fi € g* , and gAB(m) are the components of the map 
Ifi( m ) = j* ° ° j- Then 

G AB (m)» B E A = I- 1 (m)Uy), 
a point of g^. So finally the reconstruction equation may be written 

This is an equation between curves in g^. 

We will now show that this reconstruction equation above is a particular and simple case of one 
of the reconstruction equations appearing in |10j . 



22 



To do so we must introduce yet another connection, used in |10j and called there the mechanical 
connection for the G^-action. This is a connection on the principal fibre bundle M — > M/Gn, 
i.e. a G^-invariant splitting of the short exact sequence 

O^Mxg^TM^M x M/G(i T(M/G^) -> 0, 

If, as before, {Ea, E a } is a basis of g for which {Ea} is a basis for 0^, then the vector fields Xi 
together with the vector fields E a — G AB ' gAa.E>B form a basis for the set of vector fields which are 
horizontal with respect to the mechanical connection for the G^-action. We denote the latter by 
uj* 1 . Now uj^ and O are related somewhat as a connection and its vertical lift: in fact (for their 
projections onto 0„) 

n{Z v ) = i^(T(r o t)Z v ), Z v € TN^. 

We note in passing that since T(r o i)T(v) = v for any v € Nn, we can write the reconstruction 
equation as 



The reconstruction equation in |10j that we are aiming for is the third of the four, equation 
(IV. 6). It seems the one most relevant to our approach because, as Marsden et al. say, in it 
they 'take the dynamics into account', and this has been our purpose throughout. Now equation 
(IV. 6) of [10J differs from our reconstruction equation (expressed in terms of In) by having an 
additional term on the right-hand side involving the mechanical connection for the G^-action 
lo^. This arises because the authors start with a more general class of curves on M than we do. 

In order to show that our equation agrees with theirs we first show that the curve c = r o i o v n in 
M is u^-horizontal; the extra term in their equation is therefore zero in our case. By evaluating 
uj^ on the tangent to c and using the relation between uj^ and f2 we have 

c^(c) = cj^(T(t o t)u' H ) = Q(v'») = 

because v 11 is ri-horizontal. So our reconstruction equation formally agrees with equation (IV. 6) 
of Marsden et al., when we take the starting curve on M to be c: it is the particular case of that 
equation in which the curve on M is horizontal with respect to the G M mechanical connection. 

To finish the story we must also take into account the fact that equation (IV. 6) of |10j is 
presented as an equation for the reconstruction of a base integral curve of T, with momentum fM, 
from another suitable curve on M, whereas our reconstruction equation gives an integral curve of 
T on N^. But there is no real discrepancy here, because T is a second-order differential equation 
field and so knowing its base integral curves is equivalent to knowing its integral curves. Let us 
spell this out in detail. We know that if t *— > g(t) is a solution of our reconstruction equation 
then 

t » V 9 > H (t) 

is an integral curve of T in Nn. The corresponding base integral curve is 



But 



t ~ TW^ H (i))). 



r o i o ij) p = t o oi = tfj o r o l, 



so the curve t i— > ip^ t ^c(t) is a base integral curve of T. Thus the same curve in G M determines 
an integral curve of T in Nn (by its action on v H ) and the corresponding base integral curve (by 
its action on c = roioi H ] the projection of v 11 to M). 
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8 Illustrative examples 



We give two examples. In the first we derive Wong's equations using our methods. This example 
is intended to illustrate the Routhian approach in a case of some physical interest; however, we do 
not pursue the calculations as far as the consideration of the isotropy algebra and reconstruction. 
These matters are illustrated in the second example, which is more specific and more detailed, 
if somewhat more artificial. 



8.1 Wong's equations 

We discuss the generalized Routh equations for the geodesic field of a Riemannian manifold on 
which a group G acts freely and properly to the left as isometries, and where the vertical part 
of the metric (that is, its restriction to the fibres of tt m : M — > M/G) comes from a bi-invariant 
metric on G. The reduced equations in such a case are known as Wong's equations (2j 112j . 

This is of course an example of a simple mechanical system, with V = 0; we therefore adopt the 
notation of Section 7, and we will use the mechanical connection. In order to utilise conveniently 
the assumption about the vertical part of the metric g, we will need symbols for the components 
of g with respect to the invariant vector fields E a introduced in Section 2; we write 

h a b = g(E a , Eb) = A c a Afg c d- 

Since both h a b and are G-invariant functions, they pass to the quotient. In particular, the 
gij are the components with respect to the coordinate fields of a metric on M/G, the reduced 
metric; we denote by Tfj its Christoffel symbols. 

The further assumption about the vertical part of the metric has the following implications. It 
means in the first place that £g g(E a ,Eb) = (as well as C^^g(E a , Eb) = 0). Taking into account 

the bracket relations [E a ,Eb] = C c ah E c , we find that the h a b must satisfy h a dC bc + hbdC^ c = 0. 
It is implicit in our choice of an invariant basis that we are working in a local trivialization 
of M — > M/G. Then the h a b are functions on the G factor, so must be independent of the 
coordinates x l on M/G, which is to say that they must be constants. Moreover, E a , E a and A b a 
are all objects defined on the G factor, so are independent of the x\ We may write 

for some coefficients 7? which are clearly G-invariant; moreover [Xi,E a ] = jfC^Eb- We set 
liC h ac = tL; then h ac -yf b + h bc ^ b = 0. 

We are interested in the geodesic field of the Riemannian metric g. The geodesic equations may 
be derived from the Lagrangian 

L = \g a pu a vP = \g ij v i v j + \g a bv a v b = \g ij v i v j + \h ab w a w h , 

where the w a are quasi- velocities relative to the E a ; we have A b w b = v a . The momentum is 
given by p a = g a bV h = A a hb c w c , where (A a ) is the matrix inverse to (A a ). The Routhian is 

K = y ijV V - \g ah V aVb- 
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It is easy to see that X^iJZ) = gijV^ . The calculation of Xf(1Z) reduces to the calculation of 
Xi{gij) and Xi(g ab ). The first is straightforward. For the second, we note that g ab = A a Afh cd ; 
since the right-hand side is independent of the x l , so is g a b, and so equally is g ab . It follows that 

xm) = h^v k -l 1 fE c (g ab )p aPb . 

Now E c (g ab ) = —A d (g ae C de + g be C de ), from Killing's equations. Using the relation between g a b 
and h a b, and the fact that ad is a Lie algebra homomorphism, we find that 

E c (g ab ) = -A a d A b (h d ^C e cf + h^C%). 

The expression in the brackets vanishes, as follows easily from the properties of h a b- Thus the 
generalized Routh equation is 

J t (9vv j ) - f^v k = 9ij (V + l£t;V) = 

But [i a = g a bV b = A c a hb c w b ; so if we set Kfj = A b R b j, then HaRfj = hb c K£jW b . The generalized 
Routh equation is therefore equivalent to 

x* + Ti k x3x k = g™h hc K$ m x l w b . 
We also need an equation for w a : this comes from the constancy of /U a , which we may write as 

h bc j t (Alw b ) = o. 

If we are to understand this equation in the present context, we evidently need to calculate A b a . 
Now 

A b a = v*UA b a ) + v c E c (A b a ) = HtaA + v c C b cd A d a . 

It follows that 

= -hbcAiiv^fa + w^^), 

where in the last step we have again used the fact that ad is a Lie algebra homomorphism. Now 
from the skew-symmetry properties of h a b we obtain 

hcj t (A c a ) = h cd A d a (v^ b + w e C c eb ), 

and therefore 

h bc j t (A c a w b ) = h cd A d a (w c + 7 >V). 

The generalized Routh equation and the constancy of momentum together amount to the mixed 
first- and second-order equations 

■ Vj^ J x k = g m h bc K[ m x l w b 
w a + jj b x j w b = 0. 

These are Wong's equations as they are usually expressed. 
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\x 2 + \y 2 + \z 2 + \9 2 + A((sin 0)i + (cos 0)y)0. 



8.2 A Lagrangian with SE(2) as symmetry group 

We now consider the Lagrangian (of simple mechanical type) 
L 

The system is regular if A 2 ^ 1. The Euler-Lagrange equations are 

x = 0, ^r(y + (Acos9)9) = 0, -^(i + (Asm9)9) = 0, + (Asm9)z + (Acos0)y = 0, 

and the solution with (for convenience) 9$ = is 

(x(t),y(t),z(t),9(t)) 

= (x t + xq, -A sin(0 o t) + (y + A9 )t + yo,A cos(9 t) + z t + z - A, 9 t 



The system is invariant under the group SE(2), the special Euclidean group of the plane. The 
configuration manifold is R x SE(2), where x is the coordinate on R. We will use the trivial 
connection. An element of SE(2) can be represented by the matrix 



cos t) — sin t> y 
sin 9 cos 9 z 
1 



V 

The identity of the group is (y = 0, z = 0, 6 = 0) and the multiplication is given by 

(yi,zi,0i) * (y2,Z2,0 2 ) = (y 2 cos9 1 - z 2 sin(9i + y 1 ,y 2 sm6 1 + z 2 cos9i + zi,9i + 9 2 ). 
The matrices 
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form a basis for the Lie algebra, for which [ei,e2] = 0, [ei,e3] = e 2 and [e2,es] 
corresponding basis for the fundamental vector fields is 



-ex. The 



ei 



e 2 



d_ 

dy ' " z dz ' 
and for the invariant vector fields we get 
„ d .8 



8 8 8 

e 3 = -2-5- + y— + 

<9y 



dz 89 1 



e\ = cost/— — h sinf— , 



e 2 



„ 9 „ 9 

smf— + cost*—, 
oy oz 



89' 



One can easily verify that the Lagrangian is invariant. 



Before we calculate an expression for the level sets p a = fj, a , we will examine the isotropy algebra 
M of a generic point fi = /iie 1 + fi 2 e 2 + ^e 3 in g*. The relations that characterize an element 
£ = ^ei + i 2 e 2 + £ 3 e 3 of g^ are 



ev 2 = 0, eVi = 0, e x M2 - eVi = o. 
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So if we suppose that fii and fi 2 do not both vanish — we will take them from now on to be 1 
and fi respectively — then a typical element of g M is £ = ^(ei + fie 2 ). We will also set ^3 = 
for convenience. Since g^ is 1-dimensional it is of course Abelian. 

Before writing down the coordinate version of the reduced equations in the previous sections 
we made two assumptions. First, we supposed that a part of the basis of g was in fact a basis 
of g^. So from now on we will work with a new basis {E\ = e± + fie 2 ,E 2 = e 2 ,E 3 = e^}, 
with corresponding notations for the fundamental and invariant vector fields. The Lie algebra 
brackets in this basis are [E\, E 2 ] = 0, [Ei, E 3 ] = —fxEi + (1 + H 2 )E 2 and [E 2 , E 3 ] = —E\ + \iE 2 . 
The momentum vector with which we are working takes the form (1 + l-i 2 )E 1 + fiE 2 (with ^3 = 0), 
when written with respect to the new dual basis. 

The second assumption is that we use coordinates (9 a ) = (9 A ,9 a ) on G such that the fibres 
G —> G/G^ are given by 9 a = constant. Then fundamental vector fields for the G^-action 
on G are of the form K^d/d9 B . The main advantage of this assumption is that in these 
coordinates the expressions in the reduced equations became independent of the coordinates 9 A . 
This assumption is not yet satisfied in our case for the coordinates (y,z,0). The action of G^ 
on G is given by the restriction of the multiplication, i.e. by 

(yi) * (2/2, z 2 ,9 2 ) = {y 2 + yi, z 2 , 9 2 ). 

We have only one coordinate on G^, say y' . The fundamental vector fields that correspond to 
this action should be of the form Kd/dy' . However, in the new basis, vectors in g^ are of the 
form KEi, with corresponding fundamental vector fields 

So we should make a coordinate change (y, z, 9) — ► {y 1 , z' , 9'), such that 

8_ d_ d_ 

dy' dy ^ dz ' 

This can be done by putting 

y' = y, z' = z- fiy, 9' = 9. 

We will then have coordinates (y', z', 9', x, x) on iV M , and (z 1 , 9', x, x) on N^/G^. To save typing, 
we will use y and 9 for y' and 9', and only make the distinction between z and z' . 

The first goal is to solve the reduced equations on N^/G^. They are of the form 

x l = V l (xi,6 a ,xi), 
9 a = L^K^-^Af, 

For this example there is only one coordinate x on R (we are using the trivial connection on 
SE(2) x R — > R), but the coordinates 9 a on SE(2)/G fl are (z',9). The reduced second-order 
equation in x above can be derived from the Lagrangian equation in x which is simply 

x = 0. 

It is therefore not coupled to the first order equation in (z', 9), and its solution is x(t) = xot + xo. 
For the other equations, we will work first with the variables (y, z, 9), and only make the change 
to the new coordinates at the end. 
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The matrix (K^) in the above expressions is determined by the relation E a = K h a djd6 a . It is 
the lower right (2,2)-matrix of 




K 



With the trivial connection, the equations for the other variables on N^/G^ are therefore of the 
form 

z' = i 2 + (y + fiz)t 3 , 6 = l 3 . 

We can find the functions t a by solving the expressions p a = fi a for v a , with (^1,^2,^3) £ 0* of 
the form ((1 + /i 2 ), fi, 0). We get 

1 + fi 2 = y + /j,z + (A cos 9 + A/j, sin 6)6, 
H = z + (Asm6)6, 

= (A cos 6)y + (A sin 6)z + 6 - z(y + (A cos 6)6) + y(z + {A sin ( 



At t = 0, the above equations relate the integration constants and fi. We will set from now on 
yo = 1 — A60, zq = fj, and z$ = fiyo + Ayo + 6q. It is easy to see that the coordinates v a with 
respect to the basis {E{\ (with the trivial connection) are given by 

v 1 = y + z9, v 2 = z — fiy — fiz6 — yO, v 3 = 6. 

After substituting this into the equations for the level set, we obtain the expressions v a = i a as 
functions of (y,z,6). After some calculation, the reduced equations become 

A / \ 
z' = -{^(z - fiy)(sm6 - ficos6) - A(l - /j, 2 )sm6cos6 - fiA + l^A^cosd^y 

6 = -p (jxy — z + A cos 6 + Afi sin 6^j . 

Observe that we can now replace (z — /iy) everywhere by the new coordinate z', so that indeed 
the G M -coordinate y' does not appear in the reduced equations. One can verify that the solution 
of the above equations, with the integration constants determined by /i, is 

(z'(t),6(t)) = (Acos(6^t) + Atism(6 t) + (1 - A 2 )6 ,6 t) . 

We will now use the mechanical connection to reconstruct the G M -part y(t) of the solution. The 
Hessian of the Lagrangian, in the basis {X = d/dx,E a } is 

/ 1 + /i 2 /i A cos 6 + A\i sin 6 - z + fiy \ 

H 1 Asin6 + y 

A cos 6 + Afx sin 6 - z + fiy A sin 6 + y 1 - 2Az cos 6 + 2Ay sin 6 + y 2 + z 2 

\ 1 J 

The determinant of the matrix is 1 — A 2 . The vector field X° = X° = d/dx is tangent to the 
level sets and horizontal with respect to the mechanical connection £l m . 

In general, we regard the Hessian as a tensor field along the tangent bundle projection. A basis 
of vector fields along r that lie in the g-complement of is 

{^2-^2^ Es-^^iAcose + A^sme-z')^, A j . 



28 



Notice that they are all basic vector fields along r (i.e. vector fields on M). The reason is 
that the Lagrangian is of the simple type. We have seen that in that case the g-complement 
of g defines a connection on M — > M/G^. The connection tensor S7 m of the mechanical 
connection is determined by 

n m (E^) = it?, n m (K) = 0, n m (x c ) = o, n m (x v ) = o, 

where the vector fields {E&} ^at are horizontal with respect to the mechanical connection Q m 
and tangent to the level set are here 

^ = ^ ~ TT7^' E " = ^ ~ TT7 {A cos 9 + AiX sin 9 ~ z ' )& ^ 

The explicit expressions of the E% are not of direct concern, we only need to know that they 
are tangent to the level set and that they differ from E% in a vertical lift. The vertical part of 
r = i a E^ + v l Xf + T l XJ (the restriction of the dynamical vector field to N^) is here 

ft m (r) = (V + — ^i 2 + —^(Acos9 + Afisme - z')l 3 )e?- 

Before we can write down the explicit form of the reconstruction equation g~^g = f2 m (r o v H ), 
we need to find the horizontal lift v 11 of the reduced solution v = (z',9,x,x). It is the curve 
(y m , z' , 9, x, x) in whose tangent vector is horizontal with respect to the G M -mechanical 
connection. By construction this means that ^(r o v H ) should be ^-horizontal. If we write in 
general that Jj(t o v H ) = v l E\ + v 2 E2 + v 3 E%, then in order for the curve to be horizontal the 
v a must satisfy 

-y 1 = -v 2 — ^_ - v 3 (Acos9 + Ausm9 - z). 

1 + n z 1 + n z 

By expressing the v a as functions of the 9 , we find that the missing y H is a solution of 

y"H = _ A Q cos (9 t), 

from which y H (t) = — Asin(#oi) + uq. Using this y H in the reconstruction equation gives 

yi = l 1 H ^—[2 _| (Acos9 + All sin 9 - z')l 3 = 1, 

1 + \i l 1 + n l 

once we have evaluated the functions L a in terms of (y H , z', 9, x, x). So the solution through the 
identity is yi(t) = t. The y-part of the complete solution of the Euler-Lagrange equation is 
therefore 

y(t) = Vl (t) + y H (t) = -Asm(9 t) +t + y , 
as it be should for the given value of the momentum. 
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